function ejer6a
    
    A = [3 1 1; 2  6 1; 1 1  4] 
    b = [5; 9; 6]
    x = [0; 0; 0];
	er = 10^-11;
	m = 1000;
      
    fprintf('solucion de jacobi \n');
    [x1,k1] = jacobi(A, b, er, m, x)
    
    fprintf('solucion de gauss seidel \n');
    [x2,k2] = gseidel(A, b, er, m, x)
    
 endfunction   